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Abstract 

Different approaches to quantum cosmology are studied in order to deal 
with the future singularity avoidance problem. Our results show that these 
future singularities will persist but could take different forms. As an ex- 
ample we have studied the big rip which appear when one considers the 
state equation P — up with ui < — 1, showing that it does not disappear 
in modified gravity. On the other hand, it is well-known that quantum geo- 
metric effects (holonomy corrections) in loop quantum cosmology introduce 
a quadratic modification, namely proportional to p 2 , in Friedmann's equa- 
tion that replace the big rip by a non-singular bounce. However this modified 
Friedmann equation could have been obtained in an inconsistent way, what 
means that the obtained results from this equation, in particular singularity 
avoidance, would be incorrect. In fact, we will show that instead of a non- 
singular bounce, the big rip singularity would be replaced, in loop quantum 
cosmology, by other kind of singularity. 

Pacs numbers: 98.80.Qc, 04.20.Dw, 04.62,+v 

1.- Introduction — Studies of distant type la supernovae [1, 2J indicates that the 
dominant part of the energy of the universe must be gravitationally repulsive driv- 
ing our universe expanding in an accelerating way. To explain this acceleration 
one usually assumes the existence of dark energy with a negative pressure, in gen- 
eral one can assume a perfect fluid with state equation P = up, with oj < —1/3 
in order to have cosmic acceleration. Moreover, observations from WMAP indic- 
ates the value uj = —1.10 0, what means that our universe would be dominated 
by "phantom energy" (w < —1). However, the classical solutions of general re- 
lativity for a Friedmann-Robertson-Walker (FRW) model containing dark energy 
lead, in general, to future singularities |4l[5]|6l (big rip, future sudden singularities, 
etc.). Lately, a good number of papers have been dealing with the possibility of 
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avoiding these future singularities, using different approaches to quantum cosmo- 
logy like loop quantum cosmology, semiclassical gravity, modified gravity, brane 
cosmology, etc. This paper has two main objectives: The first one is to discuss 
this different approaches, to show in which way they modify the dynamics of our 
universe and to check if, effectively, they could avoid the future singularities that 
appear in classical cosmology. And the second one is to show that the modified 
Friedmann equation in loop quantum cosmology could have been obtained in an 
inconsistent way. And thus, the current statement that, in loop quantum cosmology, 
the big rip singularity is replaced by a non-singular bounce would be incorrect. 

The units used in this paper are h = c = M p = 1 being M v the reduced Planck 
mass. 



2. - Einstein cosmology — For the flat FRW spacetime filled by a perfect fluid with 
state equation P = f(p), Einstein theory, is obtained from the Lagrangian C = 
^Ra 3 — pa 3 where R = 6(H + 2H 2 ) is the scalar curvature, a is the scale factor 
and H = - is the Hubble parameter. 

This Lagrangian has been constructed in co-moving fluid coordinates (see 
and Section III C of [ 8 ]), and the energy density p has to be understood as a function 
of the scalar factor a. This relation comes from the consevation equation 

d(pa 3 ) = -Pd(a 3 ) ^p = -3H(p + P) <=> ^ = - 3 (p + P). (1) 

da a 

Since the state equation that we are studying is P = f(p) one has the differen- 
tial equation 

dp "da, (2) 



P + f(p) a 

that after integration gives p as a function of a. For example, when P = cop one 
has 

( ffl r 3(1+u) 

Pip) =Po[ — ) , (3) 

where po is the value of p when a = a^. 

The Lagrangian can be written as follows C = 3 ( d( - a ? '- — d 2 a) — pa 3 , this 



dt 

means that the same theory is obtained avoiding the total derivative, which gives 
the Lagrangian £e = —3H 2 a 3 — pa 3 . The conjugate momentum is then given by 
p = ^fi- = —6Ha 2 , and thus the Hamiltonian is 

7i E = ap- C E = -3H 2 a 3 + pa 3 . (4) 

In general relativity the Hamiltonian is constrained to be zero, which gives the 
Friedmann equation 

r2 



H = p/3, (5) 
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that together with the conservation equation are the dynamical equations that de- 
scribe the evolution of the universe. 

The Raychaudury equation is obtained from the Hamilton equation p = — ^Uf-, 
which gives 



p=-A- d -H^=-A+^ 2 > (6) 

12a 2 da 12a 2 



where we have used the conservation equation. Then from the Friedmann equation 
one easily obtains H = — |(p + P). 

A solvable example is the case of a barotropic perfect fluid with state equation 
P = top [4\, which gives 

m = sal u)t- ts p{t) = 3(1 1 u)Ht- ts r (7) 

where t s = to — 3go ( 2 1+aJ ) » being Ho = H(to) the initial condition. Then, if 
one assumes Ho > and to < — 1 one has t s > to, and thus, one has a big rip 
singularity. 

3. - Modified gravity — An alternative to Einstein cosmology is modified gravity, 
where higher-curvature terms are taken into account. The theory is based in the 
Lagrangian (see for instance EHHSHD) C>MG = f(R)a 3 - a 3 p, and to find the 
Hamiltonian formulation, one can use the Ostrogradskii's construction [12] intro- 
ducing the variables a\ = a and <i2 = a, and the momenta 

_ dC M G d 8C M G a 2fll f r>\i, _ 9C M G r 2 f /fV\ /q\ 

Pi = -9tT-J t -MT = - 6 »f ( R ) R , j* = -^- = 6a/(H). (8) 
Then, the Hamiltonian in modified gravity is given by 

UMG = Pid + P2'd - C M G = 
-6f"(R)RH + f'(R)(R - 6H 2 ) - f(R) + p)a 3 , (9) 



and the Hamiltonian constraint 7~Lmg = gives the modified Friedmann equation 

- 6f"(R)RH + f'(R)(R - 6H 2 ) - f(R) + p = 0, (10) 

that with the conservation equation p = —3H(P + p) gives the dynamics of the 
universe in modified gravity. 

Note that, the Hamilton equations di = ^f 2 , d 2 = andp 2 = 

are identities. The dynamical equation, i.e. the Euler-Lagrange equation, is p\ = 
— 8 g^ G which gives the modified Rauchaudury equation. 

Remark 0.1. A more general theory consist in modified Gauss-Bonet gravity which 
is based on the Lagrangian Cqb = f{R, G)a 3 — a 3 p, being G = 24H 2 (H + H 2 ) 
the scalar Gauss-Bonet curvature (see for instance H13V ). 
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As special case we can consider f(R) = =f — j^R 2 - Then, equation (10) 



becomes 

H 2 = - + 2a{3H 2 H + HH - -H 2 ), (11) 
3 2 

which coincides with the semiclassical Friedmann equation, obtained taking into 
account quantum effects due to a massless conformally coupled field when one 
chooses the other parameter, namely /3, equal to zero (see for instance [14J). 

For the state equation P = ujp, the case u < — 1, a > was studied in great 
detail in lTT4l . There, the main obtained result is that almost all the solution have 
future singularities in the contracting phase, having the following behavior: 

H{t) ~ - 1 p(t)~0, for t<t s (12) 

being t s the time at which the future singularity appears. This behavior means that 
if the universe is initially in the expanding phase, like nowadays, it will bounce and 
will enter in the contracting phase where it will develop the singularity described 



by (12 1. 



On the other hand, the case w < —I, a < Ois completely different, because 
now the expanding and contracting phase decouple, that is, the bounces are not 
allowed. To show this, we consider the new variable p 2 = eH where e = sign(H). 
Using this variable the modified Friedmann equation becomes 

j t (p 2 /2 + V(p, p)) = -3ep 2 p 2 + ^-(1 + u)p 

<=^. p = -d p V(p, p) - 3ep 2 p, (13) 

where V(p, p) = — §^ {p 2 + ^z)- Then, since V(0, p) = +oo this means that 
the universe cannot bounce from one phase to another. 

Once we have proved that the universe cannot bounce, the next step is to look 
for singular solution in the expanding phase and compare them with the classical 
ones (equation ([7])). To do this, we look for solutions with the following beha- 
vior at late times p(t) ~ Po{t s — t)~ u with v > 0. Inserting this solution in the 
conservation equation one obtains H(t) ~ 3(i+J)(t _t) > anc ^ fi nai ly> inserting both 



expressions in the modified Friedmann equation (equation (111) and retaining the 



leading terms, one gets v = — 4 and po = — ^^w^ > 0. Thus, 

H(t) ; r, pit) ~ . P ° , A . (14) 

Comparing this solution with (|7]) we deduce that, in that case, modified gravity 
make worse the singularities. 

To end this Section, note that for more complicated functions f(R), it seems 
impossible to perform a qualitative analysis of the dynamics. In such cases only 
numerical simulations could show the behavior of our universe at late times. 
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4.- Loop quantum cosmology — An approach to quantum cosmology that could 
avoid the big rip singularity is loop quantum cosmology, where for states which 
correspond to a macroscopic universe, such as ours, at late times the following 
effective Hamiltonian, which captures the underlying loop quantum dynamics, is 
considered ifBl [T6l [T71 

Hlqc = -ZV 8 ^^ + Vp, (15) 



where 7 = 0.2375 is the Barbero-Immirzi parameter [18] and A is a parameter 
with dimensions of length, which is determined invoking the quantum nature of 
the geometry, that is, identifying its square with the minimum eigenvalue of the 



area operator in LQG, which gives as a result A = \J ^7 (see El). In jl5l, V 
is the physical volume V = a 3 and f3 is canonically conjugate to V and satisfies 
{/3, V} = |, where {, } is the Poisson bracket, which for the canonically conjugate 
variables (a,p = -6Ha 2 ) takes the form {0, V} = f % - §g££ = -3a 2 f . 
Then, since {/?, V} = % one can conclude that j3 = "yH. 

Remark 0.2. This last statement does not seems clair in some papers of loop 
quantum cosmology. For example, in M9\l in order to define (3 the authors as- 
sert that "On the classical solution j3 is related to the scalar factor as (3 = ^H" 
and in H20\l it is stated that " The variable (3, in the limit /3 — > 0, is linked to the 
Hubble factor via the relation f3 = ^yH ". We will discuss later what really happens 
with the relation between j3 and H. 

The Hamiltonian constraint is then given by sm 2 ^f^ = £ , and the Hamiltonian 



7 

equation gives the following identity: 

V = {V,H LQC } = 
H = « p = i ar Csin (2A 7 ff)^ (16) 

Writing this last equation as follows H 2 = hl ° 2 ^f (1 — sin 2 (A/3)) and using 

the Hamiltonian constraint H lqc = h '"a^f = f one obtains the following 
modified Friedmann equation in loop quantum cosmology 



pj p c /12 pin 

being p c = ^r^- This equation with the conservation equation p = —3H(p + P) 
gives the dynamics of the universe in loop quantum cosmology. 
Here two remark are in order: 
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1. The Hamiltonian ( 15 1 can be actually constructed by using the general for- 
mulae of loop gravity that express the Hamiltonian in terms of holonomies 
hj(X) = er^^i where aj are the Pauli matrices I2lll22ll23l : 

n LQG = : , vi ^.^7r [hiWhjWhr^X) 

^ i,j,k 

x/ l T 1 (A)/ ifc (A){/ ife 1 (A),y}] + P V. (18) 

A simple calculation shows, see for instance |24ll251 . that ( 18 ) equals ( [T5] >. 

2. The old quantization of loop quantum cosmology was done using two ca- 
nonically conjugate variables, one of them was the dynamical part of the 
connection, namely c, and the other one was the dynamical part of the triad, 
namely p, (see for instance fl26ll27ll2"8l ). These variables are related with the 
scalar factor and the extrinsic curvature K = \a by the relations 

p = a 2 , c = 1 K= 'la. (19) 

Then, in order to obtain the dynamics of the universe, the following effective 
Hamiltonian was used I291I351I3T1 

-Hqlc = "^P 1/2 sin 2 (2/xc) + pp 3 / 2 , (20) 

where fi = (see l23l ) is obtained by identifying the eigenvalue ^ of 
the operator p with the minimum eigenvalue of the area operator in loop 
quantum gravity which is given by ^7. 

Using this Hamiltonian, the scalar factor satisfies the dynamical equation 

■ j nj x sin(4/xc) 

a = {a,n LCi = —p. • (21) 

2/x 7 

and, imposing once again the Hamiltonian constraint, Holc = 0, one ob- 
tains the following modified Friedmann equation 

H 2 = P -(l-^{), (22) 



3 V Pc(a) 

with critical density 



p c {a) = 5 2 2 . (23) 
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Coming back to the modified Friedmann equation in loop cosmology, (equa- 



tion ( 17 1), from the equation of the ellipsis one can see that the Hubble parameter 
belong in the interval [— p c /12, p c /12], and the energy density p in [0, p c \. Then, if 
the state equation P = f(p) is smooth enough, the functions H (t) and p(t) will be 
smooth and they also will be defined for all time, that is, there won't singularities. 
As an example, we consider the solvable case P = ujp. The dynamical equa- 



tions are now 



3 V p 
which solution is given by 



H 2 = Hl--^\, p = -3H(l + u)p, (24) 



p(t)=[^l + u) 2 (i-tf + y) , H (t) = ^(t-t)p(t) (25) 



2,/l- 



_P0 



where t = to 7^ — r^7= an d Po is the current energy density of our universe. 

V3(l+^) v /po 



Note that this solution is defined for all time and it finishes at (H = 0, p = 0), 
and its main property is that the universe remains in the expanding phase until time 
t = i. At this time it bounces and re-collapses forever and ever. 



It is clear that the behavior described by equation (25 ), is very different to the 
classical one ([7]), where the universe presents a big rip singularity. 

From this result, it seems that holonomy corrections replace the big rip singu- 
larity, which appears in classical cosmology, by a non-singular bounce. However, 
from our viewpoint, we have some objections to the way that the modified Fried- 
mann equation has been obtained. 

1. Loop quantum cosmology was built using two canonically conjugate vari- 
ables, one is the dynamical part of the connection, namely c, and the other 
one is the dynamical part of the triad, namely p, (see for instance ll26l l27l 
l28l ). These variables are related with the scalar factor and the extrinsic 
curvature K = by the relations (see Section II of ll27l and Section II 

of mn) 

p = a 2 , c = jK = ~a. (26) 

Later in 021 two new canonically conjugate variables were introduced, (V, f3) 
which are related with the standard variables through the relations V = a 3 
and f3 = -yH (see formulas 2.1 and 2.2 of [32]). Then, the loop quantum 



theory built with this two variables provides the effective Hamiltonian ( 15 1. 



However, if one starts directly from the effective Hamiltonian ( 15 ), although 
it is assumed that V = a 3 , the definition of the variable f3 comes from the 
equation V = {V^lqc}, which gives P = j\ arcsin(2A7-£f) and differs 
from the initial definition of (3. As a consequence, if one takes V = a 3 and 
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ft = J\ arcsin(2A7-ff) as a canonically conjugate variables, then the stand- 
ard variables (a, p = —6Ha 2 ) do not remain canonically conjugate, because 
now 

_ 7 / da dp da dp\ 
{a,P\ = g [opdv ~ dVd^J = 

cos(2A/3) = \A - 4l 2 \ 2 H 2 / constant. (27) 

One of the main reason against the modified Friedmann equation in loop 
quantum cosmology comes from the Legendre transformation 

U LQ c = - 2 VP-C LQC (28) 

7 

which gives, in terms of the standard variables, the following Lagrangian 

3a 3 H . . . 

L-LQC = arcsm(2A7ii) + 

7 A 

'\ - yjl - 4^X 2 H 2 ) - a 3 p, (29) 



2 7 2 A 2 

which coincides with C e for small values of H . 

It's well-known that the other current cosmological theories are built from 
two invariant, the scalar curvature R = 6 (^H + 2H 2 ^j and the Gauss-Bonnet 

curvature invariant G = 2AH 2 (^H + H 2 ^j . For example, in modified Gauss- 

Bonet gravity [13] the Lagrangian Cqb = a 3 f(R,G) — a 3 p is used, and 
semiclassical gravity, when one takes into account the quantum effects due 
to a massless conformally coupled field (see for instance [14]), is based in 
the trace anomaly T vac = aOR—^G (being a > and (3 < two renormal- 
ization coefficients). However, the Lagrangian ( [29] ) does not seem invariant, 
which is in disagreement with one of the main principles of general relativity. 

In semiclassical gravity or in the f(R, G) = § — f^R 2 theory, the modified 
Friedmann equation is given by 

H 2 = J + 2a(3H 2 H + HH - -H 2 ) - /3F 4 , 
3 2 

09 = in the f(R, G ) = *-^R 2 theory ) , (30) 

which contains higher-curvature terms like H and H. Using this equation, it 
was proved in lfl4l for the state equation P = ujp that, if — 1 < ^ < the 
universe will bounce but when it will enter in the contracting phase it will 
develop a future singularity of the form 




H(t)~— -1± Wl + i-\ 7l -, p(t)~0, for t<t s (31) 
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being t s the time at which the future singularity appears. It is clear that this 
behavior is very different to the one described in loop quantum cosmology 
(eq. (p5]>), where the universe does not develop any kind of singularity. 



On the other hand, one can adopt another different point of view. One can as- 
sume that the variables (a,p = —6Ha 2 ) are canonically conjugate, which means 
that V = a 3 and f3 = jH as we have seen at the beginning of this Section. 



Then, one has to understand the Hamiltoninan ( [15] ), not like the Hamilonian of 
the system, but as the new Hamiltonian constraint that replaces the classical one. 
From this viewpoint, taking the derivative with respect to the time of the Hamilto- 
nian constraint and finally using the conservation equation, the following modified 
Raychaudury equation will be obtained 



H ^ 



sin(2A7#) 



3sin 2 (A7iT) {3 sin 2 (\~fH)\ 



A 2 7 2 J \ A 2 7 2 J 



(32) 



For the case P = f(p) = up this equation becomes 



A__| (1 + „ )ff 5E^, (33) 

from which one deduces that H is positive, which means that H reach the value 
2^ in a finite time and thus, at that time H diverges or equivalently the scalar 

curvature R = Q(H +2H 2 ) diverges. Moreover, at that time, from the Hamiltonian 
constraint and the state equation P = up, one has p = ^2 and P = -f^. Then, 
one can conclude that, from this viewpoint, the big rip singularity is replaced by 
this other singularity characterized by a divergent scalar cuvature, but with finite 
values of the Hubble parameter, energy density and pressure. 

5.- Theory based on the reduced semiclassical Friedmann equation — This ap- 
proach was proposed by Parker and Simon in 1133113411 . its main idea is to obtain the 
derivatives of H from the classical Friedmann and conservation equations. Thus, 
once this has been done, one inserts these derivatives into the semiclassical Fried- 



mann equation obtained in modified gravity (equation ( 10). The result is a new 



modified Friedmann equation, but without derivatives on H. To be precise, we 
consider the theory f(R) = ^ — f^R 2 - Then, for the simplest case P = up we 
have H = + u)H 2 and H = |(1 + u) 2 H 3 , that once introduced in ( fll| ) 

provides 

H 2 = ^ + 9 -a{l + u)(2,u-l)H\ (34) 

If in this equation one makes the substitution H A = ^- (classical Friedmann 
equation) one will get the interesting equation 

H 2 = ^ + ^a(l + u)(3u-l)p 2 , (35) 
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which gives a way to obtain the modified Friedmann equations in loop quantum 
cosmology. Effectively, choosing a = — 4 Pc(i+^)(3^-i) one obtains the equation 



( 17 1. Then, in some sense, for the particular state equation P = cup, this approach 
could justify the modified Friedmann equation in loop quantum cosmology. How- 
ever, note that, for a more general state equation P = f(p), this method provides 



more complicated equations that (35 1, and thus, it is impossible to recover equation 

6. - Conclusions — Through this paper we have shown that it seems impossible to 
avoid completely the future singularities, which appear in phantom dark energy 
models, using different alternative approaches to quantum cosmology. We think 
the most efficient approach is semiclassical or equivalently modified gravity with 
f(R) = ^ — J2R 2 (with a > 0). In this theory the dynamics of the universe is 
drastically changed because the universe will bounce and will enter in the contract- 
ing phase where it will develop a future singularity like that described by equation 



( 12 1 or (31 ). The other important conclusion in the paper, is that the results about 
avoidance of future singularities, in particular the bounces, obtained from the mod- 
ified Friedmann equation in loop quantum cosmology (for example EU [36] 13710 
have to be revisited because this equation would not be justified. This does not 
mean that, in a more general theory where higher-curvature terms would be com- 
bined with holonomy corrections, features such as bounces may appear. But this is 
a complicated problem that deserves future investigations. 
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